sum2docx nev cpl ipr gea pn gdp ir idl ol using 描述性统计.docx,replace stats(N sd(%9.4f) mean(%9.4f) min(%9.4f) p25(%9.4f) median(%9.4f) p75(%9.4f) max(%9.4f))  title("Table 2: Summary Statistics")
global x pn gdp ir idl ol
**基准回归**
reghdfe nev cpl
est sto m1
reghdfe nev cpl,absorb(year)
est sto m2
reghdfe nev cpl,absorb(pro year)
est sto m3
reghdfe nev cpl $x ,absorb(pro year)
est sto m4
outreg2[m1 m2 m3 m4]using 基准回归.doc, append  tstat bdec(3)tdec(3)replace title("Table 3: jizhun")
**稳健性检验**
xtset pro year
reghdfe l.nev cpl $x ,absorb(pro year)
est sto m5
reghdfe nev cpl $x il cpi rd ,absorb(pro year)
est sto m6
reghdfe nev cpl $x  ,absorb(pro year) vce(robust)
est sto m7
outreg2[m5 m6 m7]using 稳健性.doc, append  tstat bdec(3)tdec(3)replace title("Table 4: wenjian")
**机制回归**
reghdfe ipr cpl $x ,absorb(pro year)
est sto m8
reghdfe gea cpl $x ,absorb(pro year)
est sto m9
outreg2[m8 m9]using 机制.doc, append  tstat bdec(3)tdec(3)replace title("Table 6: jizhi")

**内生性检验**
xtset pro year
ivreghdfe nev (cpl = l.cpl) $x , first
eststo:ivreghdfe nev (cpl = l.cpl) $x , first savefirst savefp(f)
est restore fcpl
outreg2 using 内生性分析0610.doc, append tstat bdec(3) tdec(3)
ivreghdfe nev (cpl = l.cpl) $x , first
est sto 二阶段
outreg2 [二阶段] using 内生性分析06102.doc, append  tstat bdec(3)tdec(3)


**异质性**
**高公共环保诉求响应**
egen median_pub = median(pub)
gen D_pub = pub > median_pub
label variable D_pub "高于pub中位数为1，低于为0"
reghdfe nev cpl $x  if D_pub==1,absorb(pro year)
est sto m10
reghdfe nev cpl $x  if D_pub==0,absorb(pro year)
est sto m11
**东部**
reghdfe nev cpl $x  if Pro==1,absorb(pro year)
est sto m12
reghdfe nev cpl $x  if Pro==2,absorb(pro year)
est sto m13
reghdfe nev cpl $x  if Pro==3,absorb(pro year)
est sto m14
**高研发支出**
egen median_lnrd = median(lnrd)
gen D_lnrd = lnrd > median_lnrd
label variable D_lnrd "高于lnrd中位数为1，低于为0"
reghdfe nev cpl $x  if D_lnrd==1,absorb(pro year)
est sto m15
reghdfe nev cpl $x  if D_lnrd==0,absorb(pro year)
est sto m16
outreg2[m10 m11 m12 m13 m14 m15 m16]using 异质性.doc, append  tstat bdec(3)tdec(3)replace title("Table 7: yizhi")
**经济后果检验**
reghdfe ie c.nev##c.cpl $x ,absorb(pro year)
est sto m17
reghdfe it c.nev##c.cpl $x ,absorb(pro year)
est sto m18
outreg2[m17 m18]using 经济后果检验.doc, append  tstat bdec(3)tdec(3)replace title("Table 8: jingjihouguo")


